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MATRICES 
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Abstract. We exhibit an explicit formula for the spectral density of a (large) random matrix 
which is a diagonal matrix whose spectral density converges, perturbated by the addition of a 
symmetric matrix with Gaussian entries and a given (small) limiting variance profile. 

1. Perturbation of the spectral density of a large diagonal matrix 

In this paper, we consider the spectral measure of a random matrix Df^ defined by Df^ = 
Dn + y^Xn, for Dn a deterministic diagonal matrix whose spectral measure converges and Xn 
an Hermitian or real symmetric matrix whose entries are Gaussian independent variables, with 
a limiting variance profile (such matrices are called band matrices) . We give a first order Taylor 
expansion, as e — )• 0, of the limit spectral density, as n — )• oo, of Df^. 

The proof is elementary and based on a formula given in |12j for the Cauchy transform of the 
limit spectral distribution of Df^ as n — )• oo. 



For each n, we consider an Hermitian or real symmetric random matrix Xn = [xfjjfj^i and 
a real diagonal matrix D„ = diag(a„(l), . . . , a„(n)). We suppose that: 

(a) the entries xfj of Xn are independent (up to symmetry), centered, Gaussian with vari- 
ance denoted by an{i,j), 

(b) for a certain bounded function a defined on [0, 1] x [0, 1] and a certain bounded real 
function / defined on [0, 1], we have, in the -L°° topology, 

al{[nx\,[ny\) — > cr'^{x,y) and o„([nxJ) — > f{x), 

(c) the set of discontinuities of the function a is closed and intersects a finite number of 
times any vertical line of the square [0, 1]'^. 

For e > 0, let us define, for all n. 




Date: January 19, 2013. 

2000 Mathematics Subject Classification. 15A52, 46L54. 

Key words and phrases. Random matrices, band matrices, Hilbert transform, spectral density. 
This work was partially supported by the Agence Nationale de la Recherche grant ANR-08-BLAN-0311-03. 

1 



2 FLORENT BENAYCH-GEORGES AND NATHANAEL ENRIQUEZ 

It is known, from Shlyakhtenko in [12, Th. 4.3] (see also [2j, which also provides a fluctuation 
result), that as n tends to infinity, the spectral distribution of Df^ tends to a limit fi^ with 
Cauchy transform 

C£{z) = / C£{x, z)dx, 

Jx=0 

where Cs{x, ■ ) is defined by the fact that it is analytic, maps the upper half-plane C"*" into the 
lower one C~, and satisfies the relation 

(1) Ceix, Z) = ^ . 

^ - f{x) -^Iy=o '^'^{x,y)Ce{y,z)dy 

Our goal here is to understand /Ug — /i for small values of e. Let us introduce the set T of test 
functions we shall use here. We define 

r = It I — > ^— ; z G C 
[ z — t 

Let us now define the Hilbert transform, denoted by of a function u: 

1 / N f uit) , f u(s — v) — u(s) , 

Jt&R s — t jyeR y 



Before stating our main result, let us make some assumptions on the functions a and /: 

(d) the push-forward ^ of the uniform measure on [0, 1] by the function / has a density p 
with respect to the Lebesgue measure on M, 

(e) there exists a symmetric function t( • , • ) such that for all x, y, (T^(x, y) = T{f{x), f{y)), 

(f) there exist rjo > 0,a > and C < oo such that for almost all s € M, for all t G 
[s - r/o, s + r/o], |r(s, t)p{t) - r(s, s)p{s)\ < C\t - 

Note that by hypothesis (f) and by the boundedness of the function /, the function 

s^p{s)H[T{s,-)pi-)]{s) 
is well defined and compactly supported. 

Theorem 1. Under the hypotheses (a) to (f), as e — t- 0, for all g £ T, 

g{s)dp£{s) = / g{s)dp{s)-e / g' (s) F {s)ds + o{e) , 



withF{s) ■.= -p{s)H[t{s, ■)p{-)]{s). 

As a consequence, if the function F{ ■ ) has bounded variations, then 

fig = p + edF + o{e). 

Remark 1. Roughly speaking, this theorem states that 

,. ,. spectral law (D^) — spectral law(D„) ,^ 
lim lim = di' . 

£->0 n— !>oo £ 

It would be interesting to let e and n tend to and oo together, and to find out the adequate rate 
of convergence to get a deterministic limit or non degenerated fluctuations. We are working on 
this question. 
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Remark 2. This result provides an analogue, for our random matrix model, of the following 
formula about real random variables (valid when Y is centered and independent of X): 

density^+^y(s) = densityj^ (s) + e ^ density^ (s) + o(e). 

Remark 3. In the case where X„ is a GUE or GOE matrix, the limiting spectral distribution of 
as n ^ oo is the free convolution of the limiting spectral distribution of Dn with a semi-circle 
distribution. Several papers are devoted to the study of qualitative properties (like regularity) of 
the free convolution (see [H [TJ IH 13 [6]j. Besides, it has recently been proved that type-B free 
probability theory allows to give Taylor expansions, for small values of t, of the moments of 
fit ffl for two time- depending probability measures fit and vt (see [5l[10l[9]j. Our work differs 
from the ones mentioned above by the fact that we allow to perturb Dn by any band matrix, 
but also by the fact that it is focused on the density and not on the moments, giving an explicit 
formula rather than qualitative properties. 

Proof. For all z G C^, we have 

(2) \Ce{x,z)\<^. 

Indeed, for all y, z such that z G C+, Ce{y, z) G C~. As a consequence, the imaginary part of 
the denominator of the right hand term of ([T]) is larger than '^{z). 

Hence by ([T]) and ([2]), as e — >■ 0, Ce(x, z) — > uniformly in x. 



From what precedes, 

Ce{x,z) 



z - fix) {z - fix) - e a^x, y)C,(y, z)dy)(z - fix)) 



'y 

1 



iz - /(x))2 y^^o 

1 aHx,y) 



1 

a'^ix,y)Ceiy,z)dy + o(e) 



r a^ix,y) , ^ . X 

Jy=o z - fiy) 



where each o(e) is uniform in x G [0, 1]. 

But for all a ^ b, -, X-, — rr = -, — \ ^-r ^ , I , hence since the Lebesgue measure 

' ' yz—aY^z—b) {a— by \^z—b z—a (z—a)'' J ' 

of the set {y G [0, 1] ; fiy) = fix)} is null, we have 



iz - Jy^, z - fiy) ' 7,=o ifix) - fiy))' \z - fiy) z - fix) iz - fix)) 

As a consequence, it follows by an integration in x G [0, 1] that 

where C(-) is the Cauchy transform of fi. 

Let us now recall that the push- forward of the uniform law on [0, 1] by / is the measure 
pix)dx and that a'^ix^y) can be rewritten a'^ix^y) = r(/(x), /(y)). Hence 

C,iz)-Ciz)=e[ [ {-^--^--l—^it-s)}^^^pis)pit)dtds + oie). 
JsmJtm z — t z — s (z — s) — 
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This allows us to write that for any test function g eT, 



lim^^^M^iM = A(^), 



where 



A(5) = / {9{t) - gis) - 9'{s)it - s)}^^^pis)p{t)dtds. 



'(s,t)eR2 

Note that by the Taylor-Lagrange formula, for all s, t, 



{git) - g{s) - g\s)it - s)}^^p{s)p{t) 



< Pis)pit) X ||r(-, OIIl^II/IIlo 



so that, since p is a density, by dominated convergence, 

^(5) = iToX,t)eR2-f5(^) - ^(^) - - s)}0^p{s)p{t)dsdt. 



\s-t\>V 

But by symmetry, for all 77 > 



\s-t\>V 

As a consequence, A{g) = liia^^Q A^{g), with 



Jis,t)m^^^^^^ - 9{s)}-0r^p{s)pmsdt = 0. 



^r^is) := l^^^^,.9'isf-^p{s)pit)dsdt. 



\s-t\>V 



Let us prove that almost all s G M, lim^^o / teM I^^^^M£)£W(Ji exists and that 



|s-t|>?? 



Aig) = I g\s) ( lim / r{s,t)p{s)p{t) \ 



dt. 



For 77 > and s G M, set 

^ I r{s,t)p{s)p{t) 

Set also M := \\f\\ioo. Then the support of the function p is contained in [— M, M], and so does 
the support of the function for any 77 > 0. For almost all s G [— M, M], lim^^o^jjCs) exists 
by the formula 

r(s, t)p{s)p{t) - t{s, s)pis)p{s) 



Us) = f 
Ju 



t&[s~2M,s-ri\yj[s+r],s+2M] S — t 



-dt 



and by Hypothesis (f). Moreover, for r/o as in Hypothesis (f), 

Jt=s+r) Jte[s-2M,s-r;o]U[s+r;o,s+2M] * 

< + r{s,t)p{s)p{t)dsdt 

« m JteM. 

a 770 



dt 
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Hence by dominated convergence, J^^j^ ^'(s) lim^_j.o ^J7('5)ds = limr/^o f^^^g' {s)9rj{s)ds, i.e. 

AeR U^o/*eM s-t 



\s-t\>V 



2. Examples 



□ 



2.1. Perturbation of a uniform distribution by a standard band matrix. Let us consider 
the case where f{x) = x (so that p is the uniform distribution on [0, 1]) and o"^(x, y) = l|j^_^|<^, 
where £ is a fixed parameter in [0, 1] (the width of the band). In this case, t( • , • ) = fT^( • , • ) 
and 

7a(1-s)^ 



F(.) = 1(0,1) (.) log ^ 

For small values of e and large values of n, the density pe of the eigenvalue distribution p^ of 
Z)^ is approximately 

s 1 — s J 

which means that the additive perturbation y^X^ alters the spectrum of Dn essentially by 
decreasing the amount of extreme eigenvalues. This phenomenon is illustrated by Figure [T] (where 
we ploted the cumulative distribution functions rather than the densities for visual reasons). 



d ( 

Pe{s) = p{s) + e^^(s) + o(e) = l(o,i)(s) " ^ ( 



(a) Case where n 
width I = 0.2 



4.10'' 



10"^, with 




(b) Case where n 
width t = 0.9 



4.10'' 



10" 



with 



Figure 1. Perturbation of a uniform distribution by a standard band 

matrix: plot of the functions F{ ■ ) and "^^"^ ^^^^"^ ^ (with i^Dfj( " ) and -Fd„( " ) 
the cumulative eigenvalue distribution functions of and Dn) for different val- 
ues of I. 



2.2. Perturbation of the triangular pulse distribution by a GOE matrix. Let us con- 
sider the case where p{x) = (1 — |x|)l[_i (x) and o"^ = 1 (what follows can be adapted to the 
case a'^{x,y) = i\y-x\<e, but the formulas are a bit heavy). In this case, thanks to the formula 
(9.6) of H[p{ ■ )] given~p. 509 of [llj, we get 

F{s) = (1 - |s|)l[_i,i](s) {(1 - s) log(l -s)-{l + s) log(l + s) + 2s log \s\} . 
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For small values of e and large values of n, the density pe of the eigenvalue distribution //g of 
-D^ is approximately 

d 

Pe{s) = p{s) + e—F{s) + o{e), 

which implies that the additive perturbation ^J^Xn alters the spectrum of Dn by increasing the 
amount of eigenvalues in [—1, —0.5] U [0.5, 1] and decreasing the amount of eigenvalues around 
zero. This phenomenon is illustrated by Figure [2} 




Figure 2. Perturbation of the triangular pulse distribution by a GOE 
matrix: Left: plot of the functions F{-) and ^^"^ ) ^^d„{ ) ^^j^j^ ^D^i ■ ) and 
Fd^( • ) the cumulative eigenvalue distribution functions of Df^ and Dn)- Right: 
plot of the eigenvalues histogram of and of the spectral density p of Dn- 
On the right figure, the (infinitesimal) increase of eigenvalues with respect to 
p on [—1, —0.5] U [0.5, 1] and the (infinitesimal) decrease around zero can be 
observed, in agreement with the fact that, as the left figure shows, -F' S> on 
(approximately) [—1, —0.5] U [0.5, 1] and F' ^ around zero. Both figures were 
made with the same simulation {n = 6-10^ and e = 10~^). 



2.3. Free convolution with a semi-circular distribution and complex Burger's equa- 
tion. Let us consider the case where cr^ = 1, which happens for example if the matrix Xn is 
taken in the Gaussian Orthogonal Ensemble. In this case, by the theory of free probability 
developped by Dan Voiculescu (see e.g. [13] or [H Cor 5.4.11 (ii)]), for all t > 0, 

pt = pSXt, 

where Aj is the semi-circular distribution with variance t, i.e. the distribution with support 
[—2y/t,2\/t] and density ^V4:t — x^. In this case, we know by the work of Biane [HI Cor. 2] 
that for all t > 0, pt admits a density pt- By the implicit function theorem, and the formula 
given in [51 Cor. 2], one easily sees that the function (s,t) i — > pt{s) is regular. Then, by 
Theorem [T] and the fact that the linear span of T is dense in the set of continuous functions on 
the real line with null limit at infinity, one easily recovers the following PDF, which is a kind 
of projection on the real axis of the imaginary part of complex Burger's equation given in [HI 
Intro.] 

(3) UtPtis) + l{Ptm[pti-)]is)} = 0, 

\po{s) = p{s). 

For example, if /u = Ac for a certain c > 0, then by the semi-group property of the semi-circle 
distribution [1, Ex. 5.3.26], for all t > 0, pt = Xc+t and pt(s) = 2n(c+t) V ^(^^ + ^) ~ ■s^- One can 



PERTURBATIONS OF DIAGONAL MATRICES BY BAND RANDOM MATRICES 



7 



then verify ([s]), using the formula (9.21) of H[pt{ ■ )] given p. 511 of [TT]. 

Acknowledgements. It is a pleasure to thank Guy David for his useful advices about the 
Hilbert transform. 

References 

[1] G. Anderson, A. Guionnet, O. Zeitouni. ^4?! Introduction to Random Matrices. Cambridge studies in advanced 
mathematics, 118 (2009). 

[2] G. Anderson, O. Zeitouni. A CLT for a band matrix model. Probab. Theory Related Fields 134 (2005), 
283-338. 

[3] S.T. Belinschi. A note on regularity for free convolutions, Ann. Inst. H. Poincare Probab. Statist. 42 (2006), 
no. 5, 635-648. 

[4] S.T. Belinschi. The Lebesgue decomposition of the free additive convolution of two probability distributions, 

Probab. Theory Related Fields 142 (2008), no. 1-2, 125-150. 
[5] S.T. Belinschi, D. Shlyakhtenko. Free probability of type B: analytic aspects and applications, to appear in 

American J. Math. 

[6] S.T. Belinschi, F. Benaych-Georges, A. Guionnet. Regularization by Free Additive Convolution, Square and 

Rectangular Cases. Complex Analysis and Operator Theory. Vol. 3, no. 3 (2009) 611-660. 
[7] H. Bercovici and D. Voiculescu. Regularity questions for free convolution, Nonselfadjoint operator algebras, 

operator theory, and related topics, 37-47, Oper. Theory Adv. Appl. 104, Birkhauser, Basel, 1998. 
[8] P. Biane. On the Free convolution by a semi-circular distribution. Indiana University Mathematics Journal, 

Vol. 46, (1997), 705-718. 
[9] M. Fevrier. Higher order infinitesimal freeness, arXiv. 
[10] M. Fevrier, A. Nica. Infinitesimal non-crossmg cumulants and free probability of type B, J. Funct. Anal. 258, 

2983-2023 (2010). 

[11] F. W. King. Hilbert transforms. Vol. 2. Encyclopedia of Mathematics and its Applications, 125. Cambridge 

University Press, Cambridge, 2009 
[12] D. Shlyakhtenko. Random Gaussian band matrices and freeness with amalgamation. Internat. Math. Res. 

Notices 1996, no. 20, 1013-1025. 
[13] D. Voiculescu. Limit laws for random matrices and free products. Inventiones Mathematicae, 104 (1991) 

201-220. 

Florent Benaych-Georges, LPMA, UPMC Univ Paris 6, Case courier 188, 4 Place Jussieu, 75252 
Paris Cedex 05, France and CMAP, Ecole Polytechnique, route de Saclay, 91128 Palaiseau Cedex, 
France. 

E-mail address: florent .benaychOupmc . fr 

Nathanael Enriquez, MODAL'X, 200 Avenue de la Republique, 92001, Nanterre and LPMA, Case 
courier 188, 4 Place Jussieu, 75252 Paris Cedex 05, France. 

E-mail address: nenriquez@u-parislO.fr 



